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Tunable interactions between vortices and a magnetic dipole 
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The interactions between vortices in a thin superconducting film and one magnetic dipole in the 
presence of a magnetic field applied parallel to the film surfaces are studied theoretically in the 
London limit. The dipole magnetic moment is assumed to have constant magnitude and freedom to 
rotate. The pinning potential for an arbitrary vortex configuration is calculated exactly. It is found 
that, due to the dipole freedom to rotate, the pinning potential differs significantly from that for 
a permanent dipole. In particular, its dependence on the applied field is non-trivial and allows for 
tuning of the pinning potential by the applied field. The critical current for one vortex pinned by the 
dipole is obtained numerically as a function of the applied field and found to depend strongly on the 
field. Order of magnitude changes in the critical current resulting from changes in the direction and 
magnitude of the applied field are reported, with discontinuous changes taking place in some cases. 
The effect of vortex pinning by random material defects on the critical current is investigated using 
a simple model. It is found that if random pinning is weak the critical current remains strongly 
dependent on the applied field. Possible applications to vortices pinned by arrays of magnetic dots 
are briefly considered. 

PACS numbers: 74.25.Sv, 74.25.Qt 



I. INTRODUCTION 



The study of interactions between vortices in super- 
conducting films and arrays of magnetic dots placed in 
the vicinity of the film has received a great deal of atten- 
tion lately. The magnetic, superconducting, and trans- 
port properties of a great variety of such systems have 
been reported in the literature [l], 0. The main inter- 
est in this type of system is to enhance and modify vor- 
tex pinning, and thereby increase the critical current and 
stabilize new vortex phases. From the theoretical point 
of view, the problem is to understand how the magne- 
tization of the dots influence the vortices. Calculations 
of the interactions between vortices and magnetic dots 
have been reported by several authors, using both the 
Ginzburg-Landau and London theories 0,M B 
The experimental and theoretical work carried out in the 
above cited literature deals only with dots with perma- 
nent magnetization. In this case the interactions between 
vortices and magnetic dots result from the action of the 
inhomogeneous magnetic field created by the dots in the 
superconductor. As far as vortex pinning is concerned, 
the role played by the dots is a passive one, basically to 
set up a fixed pinning potential for the vortices. A ques- 
tion that naturally arises is what happens if the dots have 
non-permanent magnetization. In this case the dots are 
expected to play an active role in vortex pinning, in the 
sense that the magnetization of each dot now depends on 
the total magnetic field acting on it, which includes the 
magnetic field generated by the vortices and the applied 
field. The interaction between vortices and dots must 
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now differ from that for dots with permanent magnetiza- 
tion, and it may be even possible to tune it by changing 
the applied field. 

This paper studies theoretically the interaction be- 
tween vortices and dots with non-permanent magnetiza- 
tion using a very simple model. Vortices in a thin super- 
conducting film interacting with one point dipole placed 
outside the film in the presence of an applied magnetic 
field parallel to the film surfaces. The dipole magnetic 
moment is assumed to have constant magnitude and free- 
dom to rotate. This model is also of experimental inter- 
est because it is feasible to fabricate magnetic dot arrays 
with freely rotating magnetic moments, as demonstrated 
recently by Cowburn, Koltsov, Adeyeye, and Welland 0. 
These authors reported on the magnetic properties of ar- 
rays of nanomagnets made of Supermalloy, each nano- 
magnet being a thin circular disk of radius R. They 
found that for R ~ 50 — lOOnm the magnetic state of 
each nanomagnet is a single domain one with the mag- 
netization parallel to the disk plane, and that the mag- 
netization can be reoriented by small applied fields. As 
pointed out in Ref. 0, each nonomagnet acts like a giant 
magnetic moment free to be oriented by magnetic fields 
acting on it. 

In this paper the pinning potential for an arbitrary vor- 
tex configuration interacting with the dipole is calculated 
exactly in the London limit, based on the solutions of the 
London equations reported in Refs.0,0. It is found that 
the pinning potential can be tuned by the applied field. 
The mechanism responsible for it is that the pinning po- 
tential depends on the dipole orientation which, in turn, 
depends on the applied field. The vortices are not influ- 
enced by the applied field because it is parallel to the film 
surfaces, and the film is thin. When a transport current 
is applied to the film, the magnetic field created by it also 
contributes to the dipole orientation. This makes the pin- 
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ning potential dependent on the transport current and, 
as shown here, has some important consequences for the 
critical current. It is also found that, in general, the pin- 
ning potential for many vortices is not simply the sum of 
the pinning potentials for each vortex. Applications for 
one an two vortices interacting with the dipole are con- 
sidered. The pinning potentials are calculated, and their 
dependence on the applied field is studied in detail. It 
is shown that the pinning potential can be tuned by the 
field over a wide range. The critical current for one vor- 
tex pinned by the dipole is then investigated. It is found 
that the critical current depends strongly on the applied 
field. Changes in the critical current by as much as one 
order of magnitude are shown to result from changes in 
the magnitude and direction of the applied field. In the 
case of a magnetic moment parallel to the film surfaces, it 
is found that discontinuous jumps in the critical current 
take place in some circumstances, as the applied field 
magnitude is changed with the orientation fixed. The 
effect of random pinning on the critical current for the 
magnetic moment parallel to the film surfaces is investi- 
gated using a simple model. It is shown that if pinning 
is sufficiently weak the dependence of the critical current 
on the applied field is essentially unchanged, except that 
discontinuities are replaced by sharp changes. The paper 
also argues that the model under consideration is rele- 
vant for vortices interacting with arrays of nanomagnets 
similar to those reported in Ref.[i|. 

This paper is organized as follows. In Sec. ^ the total 
energy for the superconductor-dipole system is obtained 
in the London limit and the pinning potential for an ar- 
bitrary vortex configuration is calculated. Applications 
to pinning of one and two vortices are also considered. In 
Sec. IIIII the critical current for one vortex pinned by the 
dipole is calculated numerically. Finally in Sec. IIVI the 
results are discussed and the conclusions of the paper are 
stated. 



II. PINNING POTENTIAL 

The superconducting film is assumed to be planar, 
with surfaces parallel to each other and to the x — y 
plane, isotropic, characterized by the penetration depth 
A, and of thickness d <C A. The film has a distribution of 
vortices, each vortex with vorticity qj = ±1 and located 
at positions rj . The magnetic dipole, m, is located above 
the film at ro = (0,0, zo > 0), has constant magnitude 
m, but is free to rotate. Here two possibilities are consid- 
ered: i) m parallel to the film surfaces ( x — y plane), and 
ii) m free to point in any direction in three-dimensional 
space. Hereafter these two possibilities are referred to as 
parallel dipole and free dipole, respectively. An uniform 
magnetic field H is applied parallel to the film surfaces. 
The film-dipole system is shown in Fig. [^,. 

In the London limit, the total energy can be written 

as 




FIG. 1: a) Superconducting film with two vortices at n, and 
r2, a magnetic dipole, m, at ro = (0,0, 20), and an applied 
magnetic field H parallel to the film surfaces, b) Film with 
a transport current J along the y-direction and one vortex. 
Transport current generates field Hj at the dipole, and force 
Fi on the vortex, both along the x-direction. H makes an 
angle a with the x-axis. 



where E sm is the energy of interaction of the dipole with 
the magnetic field created by the screening supercurrent 
induced by it in the film, E v is the total energy of the vor- 
tices in the absence of the dipole, E vm is the interaction 
energy of the dipole with the vortices, and Eh = — m • H 
is the energy of the dipole in the applied field. It is con- 
venient to use the following natural scales for physical 
quantities. Length: £ = vortex core radius. Energy: eorf, 
where eo = ( ( / ) o/47rA) 2 is the basic scale for energy/length 
of the superconductor. Magnetic moment: 4>o z o- Mag- 
netic field: <po/X 2 . 

The energy E sm is given by H U H 



(2) 



where is the field of the screening supercurrent at the 
dipole position. For d -C A 
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where mj_ is the component of m perpendicular to the 
z-direction, that is, parallel to the film surfaces. For a 
parallel dipole, E sm is a constant, that is, independent 
of the dipole orientation, since m z = 0, mi = m. For 
a free dipole, E sm depends only on the orientation of m 
with respect to the z-axis. In this case the minimum of 
E sm occurs when the dipole is oriented parallel to the 
film surfaces {m z — 0), with minimum energy E sm — 
(e dn 2 /2)(m/cj)oZo) 2 

The vortex energy E v is given by 
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where A = 2A 2 /d, and £ is the vortex core radius. The 
above expression is valid provided that | — r& |<SC A 

The vortex-dipole energy, E vm , results from the mag- 
netostatic interaction of the dipole with the magnetic 
field created by the vortices 0, El 13 > that is 



E vm = m • b , 



(5) 



in general the sum of the pinning potentials for individ- 
ual vortices. It depends on the vortex positions through 
b, which is the sum of the individual vortex fields b s , 
Eq. ©. Thus, the dipole also generates vortex- vortex 
interactions. The dependence on H obtained in Eqs. JHJ, 
and UlUfl casts it in the role of a handle that controls the 
strength and spatial dependence of U vm - The pinning 
potentials for one and two vortices are discussed next. 



where b is the field produced by the vortices at the dipole 
position. For single vortex with vorticity q located at 
position r such that r <C A, the field, denoted by b s , is 
given by 01 
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where b^ is the component of h s parallel to the film 
surfaces. For many vortices, b is given by the sum of the 
fields b s for each individual vortex. 

In summary then, the total energy can be written as 



Et = E v 



E s 



m (b + H) 
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and depends both on the vortex positions and on the 
dipole orientation. The static interaction between the 
vortices and the dipole is the total energy for the dipole 
at its equilibrium orientation. That is, for the direction 
of m that minimizes Et with the vortices held at fixed 
positions. For a parallel dipole, according to Eq. J7J), 
the equilibrium m is parallel to b^ + H, because E sm is 
independent of the dipole orientation. For a free dipole, 
E sm also contributes to the equilibrium m. In this paper 
the effect of E sm is neglected, so that the equilibrium 
m is parallel to b + H. The conditions for the validity 
of this approximation are discussed in more detail later. 
The total energy for the equilibrium m can be written as 



E eq — K + TJ 



(8) 



where U vm is the pinning potential for the vortices, given 
by 



U,, m = -m I b i 



H I +mH 
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for a parallel dipole, and 

U vm = -m[(b_L + H) 2 + (M 2 )] 1/2 + mH , (10) 

for a free dipole. The constant term mH is added to the 
definition of U vm in order that it vanishes when there are 
no vortices present, that is when b = 0. Two important 
consequences of the dipole freedom to rotate present in 
Eqs. © and l|l(J|l are the many- vortex character of U vm , 
and the non-trivial dependence of on H. The former re- 
sults because, according to Eqs. 10 and (|10|l , U vm is not 



one vortex 



1« 

vm ' 



lv 
vrn 



The pinning potential for one vortex, denoted by t/, 
follows from Eqs. J5J, and l|10[l with b replaced by b s , 
Eq. ©. Typical results for the spatial dependence of 
for characteristic values of H, assuming that H is 
parallel to the x-axis, are shown in Figs. and The 
scale for H in Ul m is the vortex field, which is bound 
by | b s z \<\ q | d/4^o x (0 O /A 2 ), and | bi |< 0.3 | q 
d/A-KZQ x (0o /A 2 ). The details of the dependence of U, 
on the spatial coordinates and on H are described next, 
i) H 3> b s . In this case, since H is much larger than the 
vortex field, the equilibrium orientation of m is paral- 
lel to H, and reduces to the pinning potential for a 
vortex interacting with a permanent dipole parallel to H. 
Assuming that H is parallel to the ^-direction, — 
—mb s x , which, according to Eqs. JBJl coincides with the 
expression obtained in Refs.JEIESIll- I n this case Ul m 
is anti-symmetric with respect to an inversion of the vor- 
tex position (r — > — r), and to the change of sign of the 
vorticity (q — * —q). For q > it has a minimum ( max- 
imum ) on the a;-axis at x = — (+)1.3zo> with minimum 
(maximum) value Uy m /eod = — (+)0.3 x A-Kq{m/ 4>qZq), 
as shown in Fig. a. 

\i)Parallel dipole. For H = 0, is given by 
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Thus, for H = 0, Ul m has circular symmetry, is the same 
for vortices (q > 0) and anti- vortices (q < 0), and is 
attractive with a repulsive core. The minimum of 
is degenerate, located on a circle of radius r = 1.3,Zo, 
centered at the dipole, and minimum value U vrn /eod = 
-0.3x47r | q | (m/<f>oZo), as shown in Fig. 0b. Fori? ^ 0, 
the spatial dependence of U], m changes smoothly with H 
between the H = and H 3> b 8 limits, as shown in Fig.0 
The minimum of Ul m occurs when the vortex {q > 0) is 
located on the negative cc-axis. In this case b^ is parallel 
to H, and — —mb s xl so that the minimum of is 
identical to that for a permanent dipole parallel to the 
positive x-direction. However, the spatial dependence of 
differs from that for a permanent dipole, as shown 
in Fig. 0c. 

\i\)Free dipole. For H = 0, U}, m is given by 
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FIG. 2: Pinning potential for one vortex, U^m (in units of eod) 
interacting with a parallel dipole. Parameters: A = 10. 0£, d = 
zo = 2. Of;, and m — 0.1(/>oZo. External field in the x-direction: 
a) Permanent dipole , b) H = , c) H = O.O20 O /A 2 . X and 
Y in units of £. 
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FIG. 3: Pinning potential for one vortex, U^m (in units of eod) 
interacting with a free dipole. Parameters: A = 10. 0^, d = 
zo — 2.0£, and m = O.l^ozo- External field in the x-direction: 
a) H = 0, b)H = O.O40 O /A 2 , c)H = O.lfo/X 2 . X and Y in 
units of £ . 



The vortex pinning potential for H = also has circular 
symmetry, and is the same for vortices and anti- vortices. 
It is purely attractive with a minimum at r = 0, and 
minimum value U^/egd — — 47r | q | (m/0o z o), a s shown 
in Fig. 0a. For H ^ 0, the spatial dependence of C/^ 
changes with H as shown in Fig. [21 The minimum of 
also occurs when is parallel to H. For q > it is 
located on the negative x-axis, at a position that depends 
on H, between x = ( H = 0) and x = — 1.3z ( -ff 3> 
6 s ). The minimum value of also depends on H, and 
decreases as H increases from U^/eod = — 47rg(m/0o-Zo) 
for H = to U^/eod = -0.3 x 47rg(m/0 o zo) for iJ > 6 s . 
The approximation that neglects the contribution of E sm 
to the equilibrium orientation of the dipole in the case of 
a free dipole is justified if b s > b' m . According to Eqs. © 
and JfJJl, this requires that r/zo <| q \ (j)QZo/m. 



B. two vortices 

The pinning potential for two vortices, U^, with vor- 
ticities q\ and q2 , located, respectively, at ri and T2, is 
given by 

U* v m = -to | b s 1± + b* 2± + H | +mff , (13) 
for the parallel dipole, and by 

U% = -to | h\ + h s 2 + H | +mtf , (14) 

for the free dipole. In Eqs. and (fT^fl bf and b| 

are, respectively, the fields of vortices qi and q2 , given by 



Eqs. @), with r replaced by ri and r 2 . In general 
is not the sum of the pinning potentials for each vortex. 
The exception occurs only for H much larger than the 
vortex fields. In this case reduces to the energy of 
interaction of two vortices with a fixed dipole, equal to 
the sum of the energies of interaction of each vortex with 
a dipole oriented along H. The physical meaning of 
can be seen more clearly by calculating the force exerted 



by the dipole on each vortex, Fj 



2 c 



Using Eqs. (T3J| and {Tl}. The result is 

db s db s 
VP. = m • 1 F m • — - 



1,2. 



(15) 



In Eq. I|15|) m e9 denotes the equilibrium magnetic mo- 
ment, that is, | m eq \— to and orientation parallel to the 
total field acting on the dipole, namely 



b s 4- h s 
, u l± + U 2_L 



H 



eq | bf + b| + H | ' 
for the parallel dipole, and 



m, 



b s + bf, + H 
b s + b| + H 



(16) 



(17) 



for the free dipole. For a permanent dipole, the expres- 
sion for Pj is identical to Eq. (|T5)l with m eq replaced by 
the permanent dipole moment. Thus, the force on each 
vortex is the same as that exerted by a permanent dipole 
with magnetic moment m eq . However, the force on one 
vortex also depends on the position of the other vortex, 
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FIG. 4: Panels a) and b): positions (in units of £) of a vortex 
pair qi = q-z — 1 pinned by the a free dipole. In a) only one 
vortex is pinned for H > O.310o/A 2 . c) magnetic moment 
orientation (s = m/m , s v = 0). Parameter values: A = 
10. 0£, d = 20 = 2.0£ , and m = </>ozo- 



through m eg . This shows that the dipole freedom to ro- 
tate generates vortex- vortex interactions. 

As an application of the above results, pinning of two 
vortices with q± = q2 = 1 by the dipole is now considered. 
The equilibrium positions of the two vortices minimize 
the total energy, given by Eq. (7J with U vm replaced by 
U^n, Eqs. (pjj and IpHjl. The energy E sm is neglected in 
what follows. 

Parallel dipole. The total energy it is minimized when 
the vortex fields h*[ ± and b| ± are parallel to each other 
and to H. For H along the x-direction, both vortices 
are located on the x axis. In this case U^, Eq. (|13fl . 
becomes = ~m(bl x + b^), which is identical to the 
energy of interaction of two vortices located on the x-axis 
with a permanent dipole parallel to the x-axis. Thus the 
equilibrium positions of the vortices are independent of 
H, and identical to those for a permanent dipole. For 
H — the minimum of the total energy also occurs when 
the two vortices are collinear with the dipole position in 
the x — y plane, but their positions are degenerate with 
respect to the orientation of the line joining the vortices. 
Free dipole. It follows from Eq. (|14fl that the equilibrium 
positions for H = are related by inversion symmetry 
with respect to the dipole, that is xi = —Xi, yi = —yi- 
This arrangement cancels the component of the vortex 
field parallel to the film surfaces, so that the dipole is 
oriented along the z-direction. When H ^ this sym- 
metry is broken, and the vortex positions depend on H. 
Numerical results are shown Fig. As H increases the 
two vortices approach the x-axis, and eventually become 
located on it (H > 0.08<?(>o/A 2 in Fig. Hb). For the pa- 



rameters used in Fig. only one vortex is pinned by the 
dipole for large H (H > 0.31^ /A 2 in Fig.g|a). The other 
vortex depins at H — O.310o/A 2 , and is repelled away 
from the dipole. The orientation of the dipole changes 
smoothly with H, being nearly parallel to the film sur- 
faces when one of the vortices depins (Fig.0]c). 



III. CRITICAL CURRENT 

Now the critical current, J c , for one vortex pinned by 
the dipole is considered. The effect of a transport current 
density, J, applied to the film is twofold. First, it exerts 
on the vortex a force = (cj)od/c)J x z. Second, it cre- 
ates a field at the dipole position Hj = (2ird/c)J x z, 
which adds to the external field H, and modifies the vor- 
tex pinning potential. The reason is that is now 
given by Eqs. © and (fTU|) with H replaced by the total 
field Ht = H 4- Hj. The critical current depends on 
the relative orientation of J and H. Here it is assumed 
that J is fixed in the positive y-direction, and that H 
points in a direction that makes an angle a with the pos- 
itive x-axis ( see Fig-^b). In this case both and H,/ 
are along the positive x-direction, and have magnitudes 
Hj = 2-ndJ/c and Fl = 4>odJ/c. In this paper J c is 
obtained by solving numerically the equations of motion 
for the vortex. It is assumed that for J = the vortex 
is pinned at the absolute minimum of U^, and that J 
increases very slowly with time. These assumptions en- 
sure that the vortex follows the position of the minimum 
of U^m — Fl x as J increases, until J reaches a value for 
which the minimum becomes unstable, and the vortex 
depins. As J increases further, the vortex velocity also 
increases. The value of J c obtained here corresponds to 
J for which the vortex velocity reaches a small value cho- 
sen for numerical convenience. The obtained J c is slightly 
larger than J for which the minimum becomes unstable. 
This is analogous to the voltage criterion in J c measure- 
ments. The values of J are, of course, limited to J < Jd, 
where Jd = c^o/(12v / 37r 2 A 2 ^) is the depairing current. 
In the results reported next, regions where J c > Jd are 
discussed for the sake of completeness. Now that Hj also 
enters in U^, it must be compared with the vortex field. 
The maximum Hj occurs for J = J Cl and can be writ- 
ten as H Jc = 0.031d/C(J c /Jd)(^o/A 2 ). For J c - J d and 
d ~ zq ~ £, Hj c is comparable to the maximum value of 
b s j_ and less than the maximum value of b s z . As discussed 
next, Hj affects the parallel dipole more than the free 
dipole. 

Typical results for the dependence of J c on a and H 
are shown in Fig. [S] and Fig. EI For H much larger than 
the vortex field and to Hj , the pinning potential reduces 
to that for a permanent dipole oriented parallel to H, 
that is, at an angle a with the x-axis. The correspond- 
ing critical currents are shown in Fig. [5] and Fig. [B] by 
the curves labeled permanent dipole. For H comparable 
to the vortex field and to Hj, the critical currents differ 
considerably from those for a permanent dipole. Most re- 
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FIG. 5: Single vortex pinned by a parallel dipole for A = 
10£, d — zo = 2£: a) and b) Critical current (J c ) vs. a; c) J c 
vs. H; d) discontinuity field Hd vs m. Labels: m in units of 
4>ozo, H in units of </>o/A 2 . 
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FIG. 6: Critical current (J c ) vs. a for one vortex pinned by 
a free dipole. Parameters: m = 0.5</>o.zo and A = 10. 0^, d — 
zo = 2.0£. Labels: H in units of <j>o/\ 2 . 

markable is J c for the parallel dipole. The J c vs. a curves 
(Figs. 0a and0b) shown sharp jumps close to a = 180°. 
The J c vs. H curves has sharp jumps and discontinuities 
(Fig. 0c). The reasons for the above described behavior 
are explained in detail next. 

Permanent dipole ( H 3> b s ,Hj). The vortex pinning 
potential for a permanent dipole oriented at an angle a 
with the x-axis has spatial dependence like that shown in 
Fig. |21 a, rotated by a with respect to the x-axis, and is 
independent of H and J. For J — 0, the vortex is pinned 
at the absolute minimum of fZj^, located at a point in 
the x — y plane defined in polar coordinates, (p,9), by 
(p = 1.3zo, 9 = a + tt). The critical current depends 
on a and to, being a linear function of m, since is 
linear in to. It is found that J c depends strongly on a, 
being largest for a — 0, and decreasing smoothly with a, 
as shown in Fig. and This results from the spatial 



dependence of U^, as can be seen for a — 0, 180°. In 
these cases the critical current can be calculated analyti- 
cally because the vortex moves only along the x direction 
as J increases. The result is J c / Jd — 4to/0o z o f° r a = 0, 
and J c / Jd — 0.4to/^>o^o for a = 180°. The origin of this 
tenfold difference can be seen in the plot of U^ n shown 
Fig. [21a. The driving force is parallel to the x-axis in 
Fig. 0a for a = 0, and antiparallel for a — 180°. As 
can be seen in Fig. 0a, the slope of the potential barrier 
is much steeper in the positive x-direction than in the 
negative one. For other values of a, the depinning pro- 
cess is more complicated because the vortex motion as J 
increases is not confined to the direction of drive. 

Parallel dipole: Typical results are shown in Fig. The 
J c vs. a curves are shown in Fig.JSJa for m = Q.25(j>oZQ, 
and in Fig.JSJb for to = O.50o^o, for characteristic values 
of H . In both cases the J c vs. a curves differ considerably 
from those for a fixed dipole, being strongly dependent 
on H . There are sharp changes in J c close to a — 180°, 
like those for m = 0.25<f> z , H = O.OO10 O /A 2 (Fig. 0a) 
and to = O.50 o z o , H = O.O10 O /A 2 (Fig. 0b). The curve 
labeled H = in Fig. 0a is the limit of J c vs. a curve 
as H — > with a fixed. The strong dependence of J c on 
H is even more evident if J c is plotted as a function of 
H for fixed a, as shown in Fig. 0c for to = 0.5(f>oZQ. In 
this case it is found that for a > 146.25° the J c vs. H 
curves have discontinuities at H = Hd- For a < 146.25°, 
the dependence of J c on H is continuous, as illustrated 
by the curves for a = 135° and a — 90°. For a = 
135°, J c undergoes a rapid change with H around H = 
0.014(/)o/A 2 , whereas for a = 90° the change in J c with H 
is much slower. It is found that the J c vs. H curves have 
no discontinuities if to is smaller than a minimum value 
which depends on a. As shown in Fig. 0d, Hd vanishes 
at the minimum to, and increases above it essentially 
linearly with to. The reasons for the sharp changes in 
the J c vs. a curves, and for the discontinuities in the J c 
vs. H curves are explained in detail in what follows. 

J c vs. a curves. The sharp jumps close to a = 180° 
result from the dependence of on J, through Hy = 
H + Hj. Basically what happens is that as J in- 
creases, Hy rotates, becoming nearly parallel to positive 
x-direction when J reaches J c . The vortex pinning poten- 
tial changes accordingly. When J reaches J c it is essen- 
tially identical to for a permanent dipole oriented 
close to the positive x-direction, giving rise to a large 
J c . To demonstrated this, the instantaneous vortex po- 
sition, which coincides with the minimum of U^ n — F^x, 
is studied as a function of J as it increases with time 
for to = 0.5 4>qZo and for typical values of a. The val- 
ues of a chosen are a = 157.5°, 135°, for which the 
J c vs. H curves have, respectively, a discontinuity at 
Hd = O.O1130o/A 2 , and no discontinuity. The trajec- 
tories in the x — y plane described by the vortex are 
shown in Figs. 0a and [TJc. The vortex x and y coor- 
dinates are shown as a function of J in Figs. 0b and 
0d. For a = 157.5° and H = O.Oll^o/A 2 < H d , and 
for a = 135° and H = 0.011^ /A 2 , O.O130 O /A 2 , J c is en- 
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hanced by almost one order of magnitude with respect 
to the permanent dipole one. In these cases the position 
of the minimum undergoes a large displacement, from 
the initial one on the right side of the dipole (A in Fig. 
01 for J = 0, to the final one, on the left side of the 
dipole (C in Fig. EJ) for J = J c . This is accompanied 
by a flip in the direction of from near the negative 
x-direction at J = to one near the positive x direc- 
tion for J = J c . For a = 157.5° and H = O.O110 O /A 2 , 
for instance, J c — l.35J d , and Hy points at 3° with the 
x-axis and has magnitude Ht = O.O740 O /A 2 . The en- 
hancement in J c results because the vortex is effectively 
pinned by a permanent dipole oriented at a small angle 
with the positive x-axis, as can be seen from the results 
for a = 157.5° and H = O.O110 O /A 2 , shown in Fig. Qb. 
Most of the displacement of the minimum takes place for 
< J/J d < 0.5. For J/J d = 0.5, H T points at 12° with 
the positive x-axis and has magnitude Ht — O.O210o/A 2 . 
The corresponding U.^ n is essentially identical to that for 
a permanent dipole oriented at 12° with the x-axis. For 
0.5 < J/Jd < 1.35, the change in H T has only a small 
effect on U^. Thus the vortex is effectively pinned by a 
permanent dipole oriented at an angle between 12° and 
3° with the x-axis. The resulting J c is comparable with 
that for a permanent dipole with this orientation. For 
a = 135° and H = O.O110 O /A 2 , O.O130 O /A 2 , a similar 
mechanism leads to the J c enhancement, as shown in Fig. 
□c andQd. For a = 157.5° and H = 0.011500/ A 2 > H d , 
and for a = 135° and H = O.O150o/A 2 , on the other 
hand, where there is no enhancement in J c , the posi- 
tion of the minimum undergoes only a small displace- 
ment, from A to B in Fig. For a — 157.5° and 
H = O.O1150 O /A 2 , it is found that the minimum becomes 
unstable at J = 0.25 J d and that the corresponding Ht 
has magnitude Ht — O.OO70o/A 2 and makes an angle of 
50° with the positive x-axis. 

J c vs. H curves. The reason for the discontinuities 
in the J c vs. H curves can be seen from the results 
for a — 157.5°. The vortex equilibrium positions for 
H = O.O1150o/A 2 > H d and for H = O.O110 O /A 2 < H d , 
shown in Figs. [7|a and [7|b, essentially coincide from A 
to B. At B the vortex depins if H = O.O1150 O /A 2 > H d , 
but if H = 0.01 l^o / A 2 < H d the vortex only depins at 
C. However, since — F^x is a continuous function, 
the positions of the minima for two H values so close to 
each other cannot differ much from one another. What 
happens is that the minimum for H = O.O110o/A 2 < H d 
becomes unstable at B, when J = 0.25J d , and a sta- 
ble minimum appears again, at a slightly larger value of 
J. This new minimum follows closely the position of the 
H = O.O110o/A 2 minimum, becoming unstable close to 
point C. However, the vortex depins when the minimum 
becomes unstable for the first time at point B. The evo- 
lution of Ul v m - F L x with J is shown in Fig. El The J = 
minimum (Fig. 0a) becomes unstable at J = 0.25 J,; (Fig. 
EJb), but for J = 0AJ d , 0.8 J d stable minima are present 
again (Fig.0c andEld). For H = O.O1150 O /A 2 > H d , the 
minimum only becomes unstable once at point C. Thus, 
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FIG. 7: Vortex pinned by a parallel dipole for m = O.50o^o, 
A = 10. 0£, and d = zo = 2.0^. a) and c): Vortex trajectories. 
Dot indicates dipole location in x — y plane, b) and d): Vortex 
coordinates x, y vs. J corresponding to trajectories shown in 
a) and c). In d) top curves represent x/£, bottom curves y/£. 
Labels: A = initial vortex position ( J = 0). B and C = 
positions where the vortex depins. H in units of 0o/A 2 
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FIG. 8: Plot of U v ^, — Flx for parallel dipole with increasing 
J. Parameters: m = O.50ozo and A = 10. 0£, d = zo = 2. Of, 
H = O.O1130o/A 2 , a = 157.5°. Labels: J in units of J d . 



the discontinuities in J c result because for H > H d the 
minimum becomes unstable twice, whereas for H < FL d 
it becomes unstable only once. For a — 135° there are no 
discontinuities in J c vs. H, because the minimum only 
becomes unstable once for all H . The same is true for 
smaller values of a. 

Free dipole: In this case the J c vs. a curves are smooth. 
The same is true for the J c vs. H curves. The results 
for m = O.250o2o, A = 10. 0£, and d — z — 2.0£ are 
shown in Fig. EJ Large changes in J c with H still take 
place, except for small a. For H = 0, J c is, of course, 
independent of a. For J c ~ J d , Ht ~ O.O30o/A 2 , and 
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the maximum vortex field is O.O80o/A 2 > Ht- As H in- 
creases, J c approaches the permanent dipole limit. Note, 
however, that the value of H needed to reach this limit is 
larger than that for a parallel dipole, due to the effect of 
the z-component of the vortex field, as discussed in Sec. 

m 



A. random pinning 



Vortex pinning by random material defects can modify 
the strong dependence of J c on H obtained above for the 
parallel dipole. The reason is that in the presence of ran- 
dom pinning the vortex equilibrium position no longer 
coincides with the minimum of — Flx. The vortex 
motion as J increases is thus changed, modifying the de- 
pendence of J c on H. It is expected that if the random 
pinning force is small compared with the dipole pinning 
force, J c remains essentially identical to that in the ab- 
sence of pinning. Here the question of how large the 
random pinning force must be in order that the strong 
dependence of J c on H is destroyed is investigated. 

Random pinning is incorporated in the calculation of 
J c by the following simple model. The random pinning 
force magnitude, F pi is assumed given, and independent 
of H and J. If the total force exerted on the vortex by the 
dipole and by the transport current, F = — VUjJ^, + Pi, 
is such that F < F p , the vortex remains pinned. If, on 
the other hand, F > F p the vortex moves, and the total 
force acting on it is F(l — F p /F). The vortex equations 
of motion are now solved including the random pinning 
force, and with the same initial conditions as above. That 
is, with the vortex pinned at the absolute minimum of 
for J — 0. The J c vs. H curves for m = O.50oz o , A = 
10.0^, d = z Q = 2.0£, and F p = 0.25eo, 0.5eo are shown in 
Fig. |5J These value of F p correspond to critical currents 
in the absence of the dipole equal to 0.16Jd, 0.32,7,2. For 
F p = 0.25e , the J c vs. H curve ( Fig. EJa ) and the J c 
vs. a curve ( Fig. |5Jc ) are similar to those for F p = 
(Figs. 0c andJSJa). The differences are that for a = 180°, 
J c is weakly dependent on H, for a = 146.25° the J c vs. 
H curve has a sharp change instead of a discontinuity, 
and for both a = 146.25°, 135° the J c vs. H curve is 
shifted to smaller H. These similarities result because 
as J increases the vortex is pinned at location close to 
the minimum of — F^x, as shown in Fig. 0b. For 
F p = 0.5eo, on the other hand, the strong dependence of 
J c on H is destroyed, as shown in Figs. 0c and|SJd. 

In summary, the strong dependence of J c on H is pre- 
served if the pinning is weak. This means that, for the 
parameters used in Fig.[5]the random pinning force must 
be such that the critical current in the absence of the 
dipole is no larger than ~ 0.2 J^. 
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FIG. 9: Vortex pinned by a parallel dipole in the presence of 
random pinning, a) J c vs. H curves for several a. b) Vortex 
trajectories with and without random pinning. Labels: A= 
vortex initial location (J = 0). B and C= locations where 
vortex depins for H and F p as indicated, c) and d) J c vs. H 
and J c vs. a curves for various pinning strengths. Parame- 
ters: m = 0.5(/>o2o, A = 10.0^, an d d — zo — 2.0£: H in units 
of <f>o I A 2 , F p in units of eo . 



IV. DISCUSSION 

The results obtained in this paper are believed to be 
representative of low-T c superconducting films with mag- 
netic dipole arrays placed on top. First, the particular 
parameters used, d ~ zq ~ £, are typical ones. For in- 
stance, in the experiments with arrays of magnetic dots 
with permanent magnetization placed on top of super- 
conducting Nb films, reported in Refill d — ^Onm ~ £. 
The magnetic dots are separated from the film by a thin 
protective layer of thickness ~ 20nm, so that the distance 
from the magnetic dipole to the film is zq ~ £. Sec- 
ond, since the vortex pinning potential depends on the 
scaled parameters m/ipozo, d/zo and H\ 2 /<j)o, many su- 
perconducting film-dipole array systems are equivalent. 
The London limit used in this paper is valid for vor- 
tices in low-T c films. However, when a magnetic dipole 
is placed close to the film, it certainly breaks down if 
the dipole field destroys superconductivity locally in the 
film. Roughly speaking, London theory is valid as long 
as the maximum dipole field at the film is less than 
the upper critical field, that is tu/zq < 4>o / (2nS : ' 2 ) , or 
m/(4>o z o) < ( z o/£) 2 /2tt- For the parameters used above 
(zq = 2£), this gives m/(0o^o) < 0.64. The values of m 
used in this paper satisfy this condition. Thus, the calcu- 
lations reported above are within the limits of validity of 
London theory. The London limit would be a better ap- 
proximation if the present calculations were carried out 
for larger values of zq/£. However, the results for J c /Jd 
would be identical to those described above if m and d 
were scaled by the same factor as zo/£- For instance, 
if zq — > 2zq, J c /Jd would remain the same if d — > 2g? 
and m — > 2m, but the upper limit of m/4>oZo for the va- 



9 



lidity of the London approximation would increase by a 
factor of 4. The present model also breaks down if m is 
sufficiently large to create vortices in the film. For the pa- 
rameters used in the present calculations, the threshold 
value of m for spontaneous vortex creation is estimated 
asm~ 0.700^0; using the results of Ref.jSg. This value 
is larger than m used here. 

The simple model discussed here is relevant to arrays 
of magnetic dots on top of thin superconducting films, 
provided that: i) the dots are sufficiently far apart to 
neglect dipole-dipole interactions between them, ii) the 
number of vortices per dot is small enough, so that each 
dot pins at most one vortex, and the vortices are far 
enough apart to neglect vortex-vortex interactions. The 
vortices present in the film must be created by an external 
field perpendicular to the film surfaces. For the parallel 
dipole, this field does not affect the dipole orientation. 
However, this is not the case for the free dipole. The 
results obtained in this paper apply to the vortices pinned 
by the free dipole if the perpendicular field used to create 
the vortices is subsequently removed. 

Now, as an example, the conditions under which the 
results described above for the parallel dipole apply to a 
system of experimental interest are examined. The sys- 
tem consists of a typical array of nanomagnets reported 
in Ref. on top of a thin superconducting film. Assum- 
ing that £ = 20nm, it follows that for d = zq = 2£, 
A = lOf ( as in Sec. HE}, d = z = 40nm, A = 200nra, 
and 4>o/\ 2 — 500G. The value m — O.50o^o follows if 
the disk radius and thickness are chosen respectively as 
R r~j 50nm and t ~ 10nm, and the disk magnetization 
is taken as M ~ 10 2 /is/(nm) 3 . If the distance between 
disks in the array is a ~ Ifim, the dipole-dipole interac- 
tion energy , Edd ~ r?i 2 /a 3 , is small compared with the 
vortex pinning potential, since Edd ~ 10 [moa^f/^J], 
with max(Uvm) = 0-3 x ^q{m/<j>QZo). The values chosen 
for the disk radius and thickness, for the magnetization, 



and for the distance between disks are typical of those 
reported in Ref. 9]. For the parameter values described 
above, the maximum vortex field at the dipole position 
is ~ 12G. Thus, the results reported in Sec. ILTTl (Fig. I5|) 
predict that for H < 12G , J c depends strongly on H, 
like in Fig. [HJc, whereas for H > 12G, J c is that for a 
permanent dipole, and depends only on a. 

In conclusion then, this paper demonstrates that the 
interactions between vortices in a thin superconducting 
film and one freely rotating dipole can be tuned by a 
magnetic field applied parallel to the film surfaces. It is 
shown that the critical current for one vortex pinned by 
the dipole can be changed from a few tenths of the depair- 
ing current to values larger than the depairing current by 
changing the applied field. For fields much larger than 
the vortex field, when the dipole moment is stuck in the 
field direction, the critical current changes continously 
by one order of magnitude when the field is rotated 180°, 
from the direction parallel to the driving force to the 
direction opposite to it. For fields comparable to the 
vortex field the critical current is very sensitive to field 
variations. For a parallel dipole, very rapid and even 
discontinuous changes by as much as one order of mag- 
nitude are obtained for the critical current in clean films. 
It is shown that weak random pinning does not modify 
substantially this dependence, except for the discontin- 
uous changes in the critical current which are replaced 
by rapid ones. It is suggested that the results apply to 
experiments on magnetic dot arrays on top of supercon- 
ducting films. 
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